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We review a recently discovered continuum limit for the one-matrix 
model which describes "causal" two-dimensional quantum gravity. The 
behaviour of the quantum geometry in this limit is different from the quan- 
tum geometry of Euclidean two-dimensional quantum gravity defined by 
taking the "standard" continuum limit of the one-matrix model. Geodesic 
distance and time scale with canonical dimensions in this new limit, con- 
trary to the situation in Euclidean two-dimensional quantum gravity. Re- 
markably, whenever we compare, the known results of (generalized) causal 
dynamical triangulations are reproduced exactly by the one-matrix model. 
We complement previous results by giving a geometrical interpretation of 
the new model in terms of a generalization of the loop equation of Euclidean 
dynamical triangulations. In addition, we discuss the time evolution of the 
quantum geometry. 

PACS numbers: 04.60.-m, 04.60.Kz, 04.60.Nc 



1. Introduction 

Two dimensional quantum gravity is an interesting playground for quan- 
tum geometry. General ideas for string theory and quantum gravity can be 
tested by exactly solvable models. Particularly for string theory the 2d 
gravity point of view has been rather fruitful. Following the seminal work 
of Polyakov et al. [T] powerful conformal field theory methods were devel- 
oped. They allowed an exact solution for a select class of observables of 2d 
gravity (see e.g. [2j), including the coupling of simple matter models. 
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Starting earlier [3J and further developed in parallel with the continuum 
methods [4J is the method of dynamical triangulation (DT). Especially (gen- 
eralized) enumeration of the DT random surfaces by matrix models proved 
fruitful, and random matrix models became an important tool in the study 
of 2d Euclidean quantum gravity coupled to certain conformal field theo- 
ries. Morever, the current understanding is that whenever the discrete and 
continuum methods can be compared the results coincide [5j. 

For a large class of observables the discrete methods have been proven 
to be more powerful. Correlators with surfaces of higher genus can be effi- 
ciently computed by matrix model techniques [6] and fixed geodesic distance 
(propertime) correlation functions can be extracted by loop equations [7j, 
transfer matrices [8j or through bijections with random trees |9j. So far 
these results have eluded the continuum methods. 

In 1998 a different theory of 2d gravity was introduced called causal dy- 
namical triangulations (CDT) [lOj . Using computer simulations the method 
has been successfully extended to 4d quantum gravity [H]. The results are 
very promising and indicate that four dimensional gravity might be nonper- 
turbatively renormalizable. The origin of the renormalizability could be a 
nontrivial fixed point scenario as described by Weinberg [l2j . 

Although similar in spirit to Euclidean DT, the continuum limit of 2d 
CDT is significantly different. The main cause that puts 2d CDT in a 
different universality class from non-critical string theory is that in CDT 
the topology of spatial slices is fixed. This makes generic triangulations 
in CDT much better behaved, since unlike in Euclidean DT the spatial 
topology fluctuations cannot dominate the continuum limit. Consequently, 
the scaling of time in Euclidean DT is non-canonical and the Hausdorff 
dimension of the quantum geometry is du = 4. The quantum geometry 
of CDT on the other hand has a canonically scaling time variable and a 
Hausdorff dimension of d/f = 2. 

Recent developments have shown that spatial topology fluctuations can 
be included while preserving the appealing features of CDT. The main idea 
is to assign a scaling coupling constant to the spatial topology change pro- 
cess |13] and that this generalized CDT model can be described by a matrix 
model [H]. 

Before coming to the matrix model we rederive the disc function of 
pure CDT, i.e. without spatial topology change, by a simple geometrical 
loop equation. Such an equation is known to be remarkably powerful. It 
allows one to compute time dependent correlators, where time is defined 
as the geodesic distance [3 [15]. At the end of this letter we derive the 
differential equation for the time dependent propagator and show that unlike 
in Euclidean DT, but typical for CDT, the scaling of time is canonical. 

After discussing pure CDT we add a term to the loop equation that 
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(a) (6) 

Fig. 1. Illustration of the two composition moves to add a triangle. The white 
dot on the right-hand-side shows the position of the mark before the triangle was 
added whereas the black dot shows the mark after the triangle was added. 



introduces spatial topology change. Only upon adding this term we can 
relate the loop equation to the Schwinger-Dyson equation of a one-matrix 
model with a linear term in the action. Remarkably, the linear term allows 
us to take a continuum limit that is very different from the well known limit 
of Euclidean DT but very similar to the continuum limit of CDT [14j. In 
fact, the continuum limit of the generalized loop equation reproduces the 
results of [13\ [TBj . Amazingly, the continuum limit of the matrix model can 
already be taken at the level of the matrix action, giving another matrix 
model that has a direct continuum interpretation |17j . 



2. Geometrical loop equations for 2D causal quantum gravity 

In this section we compute the generating function for a set M of trian- 
gulations which are similar to the original set of causal triangulations and 
which lead to the same continuum physics: Let n denote the number of 
triangles and / the number of links at the boundary (which has one marked 
link), and assume the topology is that of the disk and denote the generating 
function ^{g, x): 

oo oo oo 

Hg,x)= J2Mg,x)U= M{n,l)g^x' = Y,pi{g)xK (1) 

/,n=0 l,n=Q /=0 

The triangulations can be generated by recursively adding triangles. In 
our model there are two possible moves. Firstly, one can glue two edges of 
the additional triangle to the triangulation, one to the marked edge and the 
other one next to it in the clockwise direction (Fig. [1] (a) ) . Secondly, one can 
add a triangle by simply gluing one of its edges to the marked edge of the 
triangulation and assigning the new mark to the new edge further clockwise 
(Fig. [T](b)). Together, the two moves give the following generating equation 
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Fig. 2. Graphical representation of the loop equation ^ and 



for large n and / (see Fig. [2]), 

Mg,x)]n,i = ^[^{g,x)]n,i+gx[^{g,x)]n,i. (2) 

To keep the nice pictorial interpretation whilst making equation ([2]) exact, 
even for n,l = 1 and n,l = 0, one defines the following derivative operator, 
see e.g. [18] . 

oo oo 
1=0 1=1 

X 

The exact generating equation can now be written as 

<^>ig,x) = 1 + gx + gxdl<l>{g,x) + gx^d^<i>{g,x), (5) 

where x and dx have the clear graphical interpretation of adding and re- 
moving boundary edges. Equation ([5l) is a simple linear equation and the 
solution is given by 

gx"^ — X + g 



9,$ = -($-l), dl^ = -{^-xpi{g)-l). (4) 



where pi{g) can be determined by demanding that the singularity structure 
of ([6]) does not change discontinuously near g = 0. The poles of ([6]) are 
located at 

x±= ^^'^^~ , l/g-x-=x+. (7) 

Since the expansion oipi{g) needs to be a power series we have that pi = x_, 
hence the disc function is given by 

H9^x) = - Pi{9) = Pi{9)^ ■ (8) 

l-pi{g)x 
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Fig. 3. Illustration of the composition moves to add a double link. 



3. The continuum limit 

Using the same scaling relations as in the transfer matrix formalism of 
CDT, 
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lg-a2A/2 



-aX 



(9) 



we reproduce the continuum disc amplitude of causal quantum gravity [TU] 



WAiX) 



Wa{L) 



(10) 



4. A matrix model for generalized 2D causal quantum gravity 

To include spatial topology change we introduce a quadratic term in the 
loop equation (see Fig. E]) 

^f3{g,x) = 1 + gx + gxdl^p{g,x) + gx^d.^^p{g,x) + px^^fi{g,xf, (11) 

where /3 is a coupling constant that determines the rate of the spatial topol- 
ogy fluctuations (see Fig. . To conform with matrix model conventions it 
is useful to introduce the following notation 



wi3{g,z) 



^p{g,x = l/z) 



With these conventions the loop equation is given by 

f5wf3{g, zf - v'{z)wi3{g, z) + qpig, z) = 0, 

where 



and 



v{z) = -gz + - ^gz^i v'{z) = -g + z - gz"^ 
,z) = l-g{pi{g,l3) + z). 



(12) 

(13) 
(14) 
(15) 
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Fig. 4. Graphical representation of the loop equation pT|) . 



Written in the form of (jl3p it is easily seen that the loop equation cor- 
responds to the Schwinger-Dyson equation of a simple one-matrix model 

r /AT \ 

(16) 



Zdisc. = J Dmexp (-^[tr v{m)]j , 



where m is a Hermitian N x A^-matrix and the functional form of the po- 
tential v{m) is given by i^^- 

The solution of the loop equation (|13p is of the following well-known 
form 

Ma, z) = j^ [v'{z) - ^v'{zY-AI3qp{g,z)^ . (17) 

At this stage the solution of the disc function is still implicit since it depends 
on pi{g,f3) through (jlSp . Demanding the solution to have only one cut in 
the complex z plane gives the explicit result 

Wfsig, ^) = ^ (-5 + z-gz^ + (gz - c)^j{z-c+){z-c-)j , (18) 
where c is the solution of a third order polynomial, 

2c^-3c^ + (2g'' + l)c = g\l-2p), (19) 



and 



1 - c ± V2^{1 -c)c-g'^ 



(20) 



5. The "causal" continuum limit 



In the well known continuum limit of the one-matrix model with poly- 
nomial potential and generic coupling constants the critical value for the 
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boundary cosmological constant z coincides with the critical value of c+ 
only. As a result, this standard limit has the peculiar feature that it leaves 
a non-scaling term as a memory of the discrete theory, since 

wf'^'-^x) = wns{x) + ai W^f "-(X) + 0(a2), (21) 
where Wns{x) is a non-scaling part and 

Wt'-{X) =(^X-^ \/xWX (22) 

is the continuum disc function. Observe that the continuum disc function 
is even sub leading in the lattice cutoff a. 

In our specific model the matrix potential is such that the critical points 
of c+ and c_ coincide. This leads us naturally to the universality class 
of two-dimensional CDT implying the same scaling relations as before ([9]), 
provided one also scales the coupling constant (3 [13] : 

/? = \9sa\ c = ^e'^^. (23) 

Contrary to the standard limit of the one matrix model our "causal" con- 
tinuum limit is free from leading nonscaling contributions, 

wp{g,x) = ^WK^gXX) + 0{a^), (24) 

where W\^g^ {X) is the continuum disc function previously derived with other 
methods [T51[TB]. 

^A,,. {X) = ^^{-{X^-K) + {X-C) ^{X + CY-2gJC^ , (25) 

where C is the solution to a third order polynomial equation, 

-kC + gs = 0. (26) 

Observe that this equation is precisely the continuum limit of equation (jl9p . 
Furthermore it is interesting to note that the structure of the discrete (jlSp 
and continuum (j25p disc functions is very similar. This is not a coinci- 
dence since, as has been noticed recently, the continuum results can also be 
described by a matrix model [T7] . 

Zcont. = J DMexp (-^[tvV{M)]^ , (27) 
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with the following potential 

ViM) = AM-^M^. (28) 

In fact, it can be shown that the continuum matrix model (1270 can be 
understood as the continuum limit of the matrix model p6|) with a standard 
combinatorial interpretation [T3]. While our continuum limit Q and (|23p 
is non-standard from the matrix model point of view it is very natural from 
a CDT perspective [T3l[T6]. 



6. Time evolution 

To see that the non-critical string theory limit of the matrix model is re- 
ally very different from our "causal" continuum limit, we briefly discuss the 
so-called fixed time (geodesic distance) two- loop amplitude Gf}{li,l2; g;t). 
This amplitude is defined as the sum over all triangulations with initial 
boundary of length li and final boundary of length I2 at fixed geodesic dis- 
tance t. With this in mind, it is natural to interpret the loop equation as a 
time dependent process O ESI [19], where each addition or subtraction of a 
triangle is a "l//i-th" part of a time step. For large li we have 

1 d 

00 

+ gOpih + 1, ^2; 5; t) + 2pY.pi{g, P)Gp{h -1-2, h^g; t). (29) 

/=o 

After a "discrete Laplace transformation" this equation becomes 

d d V 1 

—Gp{z,w]g;t) = — y{-g + z - gz^ -2l3wp{g,z))Gp{z,w]g]t)\^ . (30) 

The crucial difference between this equation and similar equations in non- 
critical string theory is the scaling of time in their continuum limits. Unlike 
in non-critical string theory the continuum limit of (jSOp involves a canoni- 
cally scaling time parameter r~at, yielding 

^G^,g^ (X, y ; T) = - A [ (x^ -K+2gs W^^g^ {X)) Ga,,, (X, y ; T)] . (31) 

This is precisely the result of the propagator derived from generalized CDT 
|131 116j . Already since the inception of CDT [10] it has been known that 
the scaling of the geodesic distance is intimately related to the Hausdorff 
dimension of the quantum geometry. One can argue that the "causal" con- 
tinuum limit of the one-matrix model is better behaved since its Hausdorff 
dimension is dn = 2 instead of dn = 4 in non-critical string theory [10] . 
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7. Conclusions 

We have described a recently found "causal" continuum limit for the one- 
matrix model [H]. With the here described combinatorial interpretation 
this limit is naturally defined when inserting a linear term in the action 
with a specially chosen coefficient. We have shown here that this coefficient 
can naturally be interpreted as an additional way to add a triangle in the 
loop equations. The associated "causal" continuum limit is very different 
from the "old" double scaling limit and exactly reproduces the known results 
of causal dynamical triangulations }10lll3j . An intriguing aspect of this new 
continuum limit is that the continuum results are also described by a matrix 
model. This matrix model has of course both an expansion in the coupling 
constants gs and A, as well as a large N expansion in powers of l/N"^, which 
reorganizes the power expansions in gs and A in convergent "subseries". 
By comparing with the generalized causal dynamical triangulation model 
the powers of jV"^''"'"^ in the large N expansion can be identified with the 
continuum causal dynamical surfaces of genus h [T71 [Hj . In this sense our 
"causal" continuum limit leads to a picture that is much closer in spirit to 
the original idea by 't Hooft [20] for QCD. 

As a complement to previous results we also derived the differential 
equation for the fixed geodesic distance two-loop amplitude. The observed 
canonical scaling of time is directly related to the fact that the quantum 
geometry obtained from our new "causal" continuum limit has Hausdorff 
dimension dn = 2 instead of dn = 4 for the "old" continuum limit of the 
one-matrix model. 

Until now, only the transfer matrix formalism has been available for 
analytical computations in causal dynamical triangulations. Here we ex- 
tended the available tools and presented new, more powerful, matrix model 
and loop equation methods. Importantly, these new methods allow one to 
analytically study simple matter models coupled to two-dimensional causal 
quantum gravity. Several of these models are currently under investigation. 
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